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Determination of a s (M^): a conformal mapping approach 
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We discuss a new class of expansions in perturbative QCD, based on the technique of conformal mappings of 
the Borel plane, and apply them for the determination of a s from the hadronic decays of the r lepton. Using the 
expansion up to fifth order in the MS scheme, the method leads to the prediction a s (M^) = 0.320 ± 0.011. 



1. INTRODUCTION 

The determination of the strong coupling a s 
is one of the most important tests of QCD. The 
recent determinations at various scales are in an 
impressive agreement among each other, leading 
to the world average PQ 



a s {M 2 z ) = 0.1184 ±0.0007. 



(1) 



The hadronic decays of the r lepton provide one 
of the most precise ways of extracting a s . The 
recent calculation of the Adler function to four 
loops L 2 triggered several new determinations of 
a s (M 2 ). One may note that the average reported 
in[l] 



a s (M T 2 ) = 0.330 ±0.014, 



(2) 



leads to the value a s (M§) = 0.1197 ± 0.0016, 
slightly higher than the global average ([lj. 

The discrepancy between the so-called contour- 
improved perturbation theory (CIPT) [3] and the 
more usual fixed-order expansions (FOPT) ap- 
pears to be the largest systematic theoretical un- 
certainty in the determination of a s {M 2 ). The 
recent works show that this discrepancy 

did not go away by adding the presently known 
higher-order terms in the expansion of the Adler 
function. An attempt to understand the differ- 
ence between CIPT and FOPT was performed in 
[7], based on previous studies [5], which exploit 
the information about the high-order behaviour 



of the perturbation expansion, encoded in the sin- 
gularities of the Adler function in the Borel plane. 
Essentially, the usual series in powers of the Borel 
variable is replaced by an improved expansion in 
powers of an "optimal" variable that maps the 
Borel cut-plane onto a disc. In the present con- 
tribution we extend the investigation by consid- 
ering other conformal mappings, useful when the 
known nature of the dominant singularities of the 
Borel tranform is also incorporated. We illustrate 
the usefulness of the new expansions using some 
realistic models of the Adler function, and apply 
the method for a new determination of a s (M 2 ). 

2. ADLER FUNCTION 

The determination of a s (M 2 ) is based on the 
the evaluation of the integral 



2wi 



ds ~ 
-oj(s)D(s), 
s 



(3) 



where U)(s) = 1 - 2s/M 2 + 2{s/M^f - (s/M 2 ) 4 
and the Adler function D(s) is expressed in per- 
turbative QCD as 



D(s) 



n>l 



(4) 



* Speaker 



where a s (/i 2 ) = a s (zx 2 )/7r. The first coefficients 
K n calculated in the MS scheme are K\ = 1, 
K 2 = 1.64, K 3 = 6.37 and K 4 = 49.08. 
For the next term the choices K 5 — 283 and 
K$ = 275 were made recently in [5] and [B|. Fi- 
nally, K n {— sf/j 2 ) depend on the coefficients of the 
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renormalization-group (RG), and the renormal- 
ization scale /i 2 is chosen as /i 2 = — s in CIPT 
and = M 2 in FOPT. 

Expansions alternative to (@| were proposed in 
the literature from various motivations. In the 
present work, the objective is to include theoret- 
ical knowledge about the high order behaviour of 
the series. From particular classes of Feynman 
diagrams it is known that K n ~ n\, so the renor- 
malized perturbation series (U) is divergent and is 
usually assumed to be an asymptotic series. From 
independent arguments it is known that correla- 
tion functions like D, regarded as functions of a s , 
are singular at a s — 0. For QED these facts are 
known since a long time [5], but they do not af- 
fect the phenomenological predictions since the 
coupling is very small. By contrast, for a large 
coupling like a s (Af 2 ) in QCD the consequences 
are nontrivial. 

The information about the high-order be- 
haviour of the series ((U) is included in the singu- 
larities of the Borel transform B(u), defined by 
the series 

B(u) = J2bnU n , b n = ^±\, n>0, (5) 

where /3o is the first coefficient of the RG f3- 
function. According to the present knowledge, 
B(u) has singularities on the real axis for u < — 1 
and u > 2, known as ultraviolet (UV) and in- 
frared (IR) renormalons, respectively. 

The expansion Q can be formally obtained 
from B(u) by means of an integral of Borel- 
Laplace type. The recovery of the function D{s) 
is actually ambiguous: there are many integral 
representations admitting ^ as an asymptotic 
expansion (for a recent discussion see [TO])- As 
shown in [11], the principal value prescription 

oo 

D(s) = — PV / e"^W B(u)du (6) 
Po J 



yields a function D(s) satisfying to a large ex- 
tent the general analyticity requirements in the 
s-plane, and we shall adopt here this definition. 



3. CONVERGENCE ACCELERATION 
BY CONFORMAL MAPPINGS 

Due to the first UV renormalon at u = — 1, the 
series ([5]) converges only in the disc \u\ < 1. The 
domain of convergence of an expansion and its 
convergence rate can be increased by expanding 
the function in powers of a different variable. The 
following two lemmas |12|13] show that one can 
find an optimal variable by means a conformal 
mapping of the u-planc. 

Lemma 1: Let T>\ and T> 2 be two simply- 
connected domains in the complex u-plane, with 
T>2 C T>\. Consider the conformal mappings 
Z\ = zi(u) : T>i — s- V = {zi : |zi| < 1} and 
z 2 = z 2 {ti) ■ T>2 — > V = {z 2 ■ \z 2 \ < 1} such that 
zi(0) = and z 2 (0) = 0. Then \zi(u)\ < \z 2 (u)\ 
for all u G T> 2 . 

Lemma 2: Let T>\ and T> 2 be the domains defined 
in Lemma 1, and B(u) holomorphic in T>\. De- 
fine the expansions B(u) = X)o° Cn,i(zi(u)) n and 
— X)o° c n,2{z 2 {u)) n , convergent in the unit 
discs \z\\ < 1 and \z 2 \ < 1, respectively. Then, at 
large n one has 

^s^r < i. (?) 

From this inequality it follows that the best 
asymptotic convergence rate is obtained with the 
variable that maps the whole holomorphy domain 
onto the unit disc. For the Adler function, assum- 
ing that there are no other singularities except the 
cuts along the real axis for u < — 1 and u > 2, 
the "optimal" conformal mapping is [8] 

yT+H- y/l-u/2 

W(Uj = ; :. (8) 

% /TTH+ y/l-u/2 

It follows that the expansion 

B(u) = J2 d nW n , (9) 

n>0 

where w = w(u), converges in \w\ < 1, i.e. in the 
whole cut u-plane, and has the best asymptotic 
convergence rate, compared to other conformal 
mappings which map only a part of the cut plane 
onto the unit disc. This led us to define the new 
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Table 1 

The quantity 5^ for the model defined in 5], calculated for a s (M^) — 0.34 with the standard and 
modified CI and FO expansions truncated at the order N. Exact value = 0.2371. 



TV 


CIPT 


FOPT 


CI W12 


FO w 12 


CI U>i 3 


FO w 13 


CI Wloo 


FO wioo 


CI W 2 3 


FO W 2 3 


2 


0.1776 


0.1692 


0.1977 


0. 


.2228 


0.2070 


0. 


.2203 


0. 


.1883 


0. 


.2524 


0.2123 


0.2099 


3 


0.1898 


0.2026 


0.2009 


0. 


.2460 


0.2030 


0. 


.2440 


0. 


.1975 


0. 


.2530 


0.2028 


0.2437 


4 


0.1983 


0.2200 


0.2263 





.2463 


0.2194 


0. 


.2460 


0. 


.2288 


0. 


.2465 


0.2206 


0.2463 


5 


0.2022 


0.2288 


0.2290 


0. 


.2440 


0.2268 


0. 


.2423 


0. 


.2310 


0. 


.2427 


0.2292 


0.2423 


6 


0.2046 


0.2328 


0.2324 


0. 


.2484 


0.2306 


0. 


.2421 


0. 


.2321 


0. 


.2431 


0.2319 


0.2449 


7 


0.2046 


0.2342 


0.2339 


0. 


.2536 


0.2331 


0. 


.2457 


0. 


.2333 


0. 


.2454 


0.2345 


0.2502 


8 


0.2017 


0.2353 


0.2339 





.2505 


0.2343 


0. 


.2484 


0. 


.2341 


0. 


.2471 


0.2347 


0.2476 


9 


0.2004 


0.2367 


0.2341 


0. 


.2431 


0.2348 


0. 


.2457 


0. 


.2346 


0. 


.2465 


0.2347 


0.2377 


10 


0.1842 


0.2390 


0.2351 


0. 


.2420 


0.2348 


0. 


.2394 


0. 


.2348 


0. 


.2436 


0.2353 


0.2337 


11 


0.1962 


0.2402 


0.2359 


0. 


.2406 


0.2348 


0. 


.2352 


0. 


.2349 


0. 


.2399 


0.2348 


0.2335 


12 


0.1123 


0.2436 


0.2362 


0. 


.2298 


0.2351 


0. 


.2349 


0. 


.2349 


0. 


.2370 


0.2374 


0.2262 


13 


0.2629 


0.2408 


0.2362 


0. 


.2229 


0.2355 


0. 


.2341 


0. 


.2349 


0. 


.2356 


0.2348 


0.2226 


14 


-0.2915 


0.2575 


0.2364 


0. 


.2242 


0.2361 


0. 


2303 


0. 


.2349 


0. 


.2354 


0.2395 


0.2314 


15 


1.1011 


0.2170 


0.2367 


0. 


.2173 


0.2366 


0. 


.2277 


0. 


.2350 


0. 


.2357 


0.2356 


0.2365 


16 


-3.362 


0.3818 


0.2368 


0. 


.2102 


0.2369 


0. 


.2305 


0. 


.2351 


0. 


.2360 


0.2343 


0.2374 


17 


9.5931 


-0.1881 


0.2368 


0, 


.2176 


0.2372 


0. 


.2356 


0. 


.2352 


0. 


.2360 


0.2533 


0.2512 


18 


-31.52 


2.144 


0.2368 


0. 


.2201 


0.2373 


0. 


.2371 


0. 


.2354 


0. 


.2359 


0.1926 


0.2665 



perturbative expansion [8] 



n>0 



Wn(s) = -3-PV 

Po 



e (w(u)) n du. 



(10) 



(11) 



By construction, the series (|10[) , when reexpanded 
in powers of a s , reproduces the coefficients K n 
known from Feynman diagrams. Moreover, under 
certain condition, the expansion (|10[) converges in 
a domain of the s-plane [8]. On the other hand, 
as shown in [5], the expansion functions W n are 
singular at a s — 0, resembling the expanded func- 
tion D itself. 

In the particular case of the Adler function in 
massless QCD, the nature of the first singulari- 
ties of the Borel transform is also known: near 
u = —1, B(u) ~ (1 + u) 11 and near u = 2, 
B(u) ru (1 — u/2)~ 12 . The exponents 71 and 
72 are positive and are given in |5|14|15) . As 
explained in 16 8] , this information can be ex- 
ploited to further improve the convergence, by ex- 
panding the product of B(u) with suitable factors 



that "soften" the dominant singularities. The ef- 
fect of a mild singularity, which vanishes instead 
of exploding at u = — 1 or u = 2, appears only 
at large orders in an expansion. So, at low or- 
ders we can expand in powers of variables that 
account only for the next branch-points of B(u). 
In general, we consider the functional 



w jk (u) = 



y/1 + u/j - y/1 - u/k 
^l + u/j + ^l-u/k' 



(12) 



which map the it-plane cut along u < —j and 
u> k onto the disc \wjk \ < 1 in the plane Wjk = 
Wjk(u). We take j > 1 and k > 2, the lowest 
values leading to the optimal mapping ([5]). 

As discussed in [7], in contrast with the op- 
timal conformal mapping, the "singularity soft- 
ening" procedure is not unique. Numerically, it 
turns out to be convenient to define the expansion 



S jk (u)B(u) = J2< k (tijk(u)) n , 



(13) 



n>0 



2 The mapping wiooi suggested in |17| , was applied in de- 
tail in |18) , and the mapping WJ13 was considered in |19| . 
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Figure 1. Real part of the Adler function of the model [5], calculated along the circle s = M% exp(i</>) 
with the perturbative expansions for N = 5. Left panel: CI expansions. Right panel: FO expansions. 
The exact function is represented by the solid line. For details see |7|13j . 
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Figure 2. Real part of the Adler function calculated along the circle s — M% exp(i<f)) with the perturbative 
expansions for N = 18. Left panel: CI expansions. Right panel: FO expansions. The exact function is 
represented by the solid line. The standard CIPT and FOPT exhibit big oscillations and are not shown. 



where the compensating factor Sjk{u) is a simple 
expression vanishing at u = — 1 and u = 2: 




The exponents j[ = 71 (1 + Sji) and 7;, = 
72(1 + ^/02)1 where 5^ is Kronecker's function, are 
taken such as to reproduce the nature of the first 
branch-points of B(u). 

By combining the expansion ()13[) with the def- 
inition ([6]), we are led to the class of expansions 

5(.s)-^4 fc Wf( S ), (15) 

n 



Wi k (s) - ^PV / e-ra du . (is) 

Po J bj k (u) 


Strictly speaking, for a fixed pair (J, k) the ex- 
pansion (|13|) converges only in the disc \wjk\ < 
mm[\wjk(— 1)|, |tiJjfc(2)|]. In particular, for j = 1 
the expansions diverge for \u\ > 2 if A; > 2 (for 
k = 2 the expansion converges outside the disc 
|u| = 2, but not on the real axis j8]). However, the 
expansion (TT5)) enters the Laplace-Borel integral 
© where, especially for small couplings a s (s), the 
contribution of high values of u is suppressed, re- 
ducing the effect of the lack of convergence. On 
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Table 2 



The quantity 5^ for a modified model with g?2> r = lj calculated for a s (M^) = 0.34 with the standard 
and modified CI and FO expansions truncated at the order N. The rows for N < 5 are identical to those 
in Table [TJ Exact value 5^ = 0.2102. 



N 


CIPT 


FOPT 


CI Wi2 


FO w 12 


CI U>i 3 


FO w 13 


CI Wloo 


FO Wiao 


CI U>23 


FO 1^23 


6 


0.2041 


0.2318 





.2263 


0. 


.2493 





.2271 


0. 


.2420 


0. 


.2284 


0. 


.2431 





.2260 





.2454 


7 


0.2041 


0.2290 





.2201 


0. 


.2628 





.2220 


0. 


.2481 


0. 


.2230 


0. 


.2472 





.2174 





.2580 


8 


0.2023 


0.2213 





.2202 


0. 


.2756 





.2164 


0. 


.2595 


0. 


.2182 


0. 


.2541 





.2136 





.2734 


9 


0.2037 


0.2110 





.2175 


0. 


.2742 





.2143 


0. 


.2686 


0. 


.2154 


0. 


.2608 





.2138 


0. 


.2706 


10 


0.1924 


0.2032 





.2055 





.2709 


0. 


.2144 


0. 


.2651 


0. 


.2146 


0. 


.2629 





.2115 





.2517 


11 


0.2124 


0.2004 





.1982 


0. 


.2905 





.2136 


0. 


.2504 


0. 


.2146 


0. 


.2578 





.2068 





.2531 


12 


0.1412 


0.2071 





.2007 


0. 


.3063 





.2111 


0. 


.2406 


0. 


.2148 


0. 


.2468 





.2081 





.2627 


13 


0.3121 


0.2117 





.2022 





.2820 


0. 


.2086 


0. 


.2449 


0. 


.2149 


0. 


.2340 





.2060 





.2133 


14 


-0.2105 


0.2344 





.2001 


0. 


.2666 





.2074 


0. 


.2459 


0. 


.2146 


0. 


.2239 





.2124 


0. 


.1338 


15 


1.2336 


0.1934 





.2009 


0. 


.2865 





.2079 


0. 


.2176 


0. 


.2142 


0. 


.2187 





.2087 





.1192 


16 


-3.147 


0.3500 





.2044 


0. 


.2562 





.2091 


0. 


.1676 


0. 


.2136 


0. 


.2175 





.2073 





.0930 


17 


9.948 


-0.2333 





.2059 


0. 


.1822 





.2102 


0. 


.1355 


0. 


.2130 


0. 


.2175 





.2275 


-0. 


.0415 


18 


-30.94 


2.084 





.2058 


0. 


.1722 





.2107 


0. 


.1345 


0. 


.2124 


0. 


.2159 





.1617 


-0. 


.1019 



the other hand, if j > 1 the expansion (fl~3|) does 
not converge for u > 1, and the effect of this di- 
vergence is expected to become visible earlier. By 
inserting into ([TBI the coupling a s (s) calculated 
by solving the renormalization group equation for 
s along the circle defined in the integral ([3]), we 
obtain the " countour-improved" (CI) version of 
the new expansions. The "fixed-order" (FO) ver- 
sion can be obtained in a straightforward way [7] . 
The expansion functions W* 2 were investigated 
in [7]. In what follows we shall investigate also 
the expansion in terms of W^ 3 , W^°° and W 23 . 

4. MODELS 

For testing the convergence of the various ex- 
pansions we consider a class of models of the type 
proposed in [5] , which parametrize the Borel func- 
tion B(u) in terms of a few UV and IR singular- 
ities and recover the Adler function by means of 
© . The free parameters are fixed by reproducing 
the known values of the coefficients K n for n < 4. 
The specific model proposed in [5] (which uses 
as input also K 5 = 283), leads to a rather large 
residue, d 2 R = 3.13, for the first IR renormalon. 
To avoid any bias related to this large contribu- 
tion, we investigated also alternative models hav- 



ing a smaller residue, <i 2 R = 1 (details and more 
results will be given in [T5]). 

For illustration we give in Table [1] the values 
of the integral 5^ defined in (|3|) for the model 
[5] , calculated with the standard and the modified 
CI and FO expansions, as a function of the per- 
turbative order TV. As discussed in [5] the stan- 
dard CIPT gives values systematically lower than 
the true result, while the standard FOPT gives a 
better approximation. The new CI expansions 
approach the exact value for increasing N (devi- 
ations appear only for the mapping u>23, for the 
reasons discussed above). The new FO expan- 
sions give in general a less accurate description, 
but with the mappings W13 and wioo convergent 
approximations are obtained also in the FO case. 

The numbers in Table [1] can be understood 
from Figs. 1 and 2, which show the real part 
of the Adler function along the circle \s\ = 
M 2 exp(i</>), calculated with the CI and FO ex- 
pansions truncated at N = 5 and N = 18, re- 
spectively. The CI expansions give a precise ap- 
proximation of D(s) along the whole circle (only 
the mapping W23 shows signs of divergence for 
N = 18, as expected). In the FO case the de- 
scription gradually deteriorates near the timelike 
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axis (4> = 0) due to the poor convergence of the 
expansion of ot s (s) in powers of a s (M^) in this 
region [3J. Note that the coupling a s (s) is cal- 
culated in the same way in the exact representa- 
tion ([6]) and the CI expansion functions (fT6]l . as 
the exact solution of the RG equation in terms 
of a s (M%), while in the FO version a s (s) is ex- 
panded in powers of a a (M%). Of course, if the CI 
expansions become very precise at high N, the 
FO expansions are expected to converge too (this 
is indeed visible in the right panel of Fig. 2 for 
the mappings W13 and u> loo ). 

In Tabled we present similar results for an al- 
ternative model. By construction the first five 
rows in Tables Q] and [D are the same. The CI 
expansions based on the mappings u>i2, u>i3 and 
Wioo approach at large N the exact value also in 
this case. In the FO case the description is less 
precise and, for the values of N considered, only 
the expansion based on wioo exhibit convergence. 

Table 3 

Expansion functions a s (M%) 

0.3198 (0.3200) 
W* 3 0.3212 (0.3212) 

Wl°° 0.3187 (0.3187) 

0.3197 (0.3198) 



5. DETERMINATION of a s 

We apply now the new CI expansions for a de- 
termination of a s (M^). Note that we do not use 
the models discussed in the previous Section, but 
only the first four K„ known in the MS scheme [2] , 
and an estimate of if 5. The results are given in 
Table [3J The first numbers are obtained with the 
phenomenological value 5^ = 0.2402 and K5 — 
283 [5], those in parantheses with 6^ = 0.2038 
and K 5 = 275 [5] . With the same input the stan- 
dard CIPT gives a s (M?) = 0.3425 (0.3421). 

The values in Table [3] are very close to each 
other, so the uncertainty related to the choice of 
the expansion is very small. Taking the average 
and adding the error according to the analysis in 
leads to 

a s (M?) =0.320 ±0.011. (17) 



The average, consistent with the determination 
based on standard FOPT [5] , is smaller by about 

0. 02 than the predictions [416] of the standard 
CIPT. Sorting out this discrepancy is important 
for the precision tests of QCD in r decays. 
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